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Abstract We have carried out a wide study of shadowing and antishadowing effects on J/y production in 
dAu collisions at ^^SNN = 200 GeV. We have studied the effects of three different gluon nPDF sets, using 
the exact kinematics for a 2 — > 2 process, namely g + g^J/xff + g as expected from LO pQCD. We have 
computed the rapidity dependence of /?cp and R^au for the different centrality classes of the PHENIX data. 
For mid rapidities, we have also computed the transverse-momentum dependence of the nuclear modification 
factor, which cannot be predicted with the usual 2 — > 1 simplified kinematics. All these observables have been 
compared to the PHENIX data in dAu collisions. 
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1 Introduction 

At high temperature and densities, QCD predicts the existence of a deconfined state of matter, the Quark- 
Gluon Plasma (QGP) which is expected to be produced in relativistic nucleus-nucleus (AB) collisions. For 
thirty years, charmonium production in hadron collisions has been a major subject of investigations, on both 
experimental and theoretical sides. J /iff production should indeed be sensitive to the QGP formation, by a 
process analogous to Debye screening of electromagnetic field in a plasma [1]. A significant suppression of the 
J/xff yield was observed at SPS energy by the NA50 experiment [2], and at RHIC by the PHENIX experiment 
in AuAu [3] and CuCu [4] collisions at = 200 GeV. In 2010 and 201 1, data have been taken at the LHC 

in PbPb collisions at ^/snn = 2.76 TeV, where the // Xjf has also been found to be suppressed [5, 6, 7]. 

However, the interpretation of the results obtained in AB collisions relies on a good understanding and a 
proper subtraction of the Cold Nuclear Matter (CNM) effects which are known to already impact the 7/1// 
production in proton (deuteron)-nucleus {pA or dA) collisions, where the deconfinement cannot be reached. 
Experiments on dAu collisions at RHIC [8, 9] have indeed revealed that CNM effects play an essential role 
at y^^Aw = 200 GeV in the production of J/xjf as well as of Y (see e.g. [10]). In particular, the shadowing of 
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the initial parton distributions due to the nuclear environment and the nuclear absorption resulting from the 
breakup of the cc pair by its multiple scattering with the remnants of the incident nuclei have a significant 
impact. 

Previous studies [11, 12, 13] have also shown that the //i// partonic -production mechanism affects the 
way to compute the nuclear shadowing and thus its expected impact on the J /^f production. Most studies on 
the 7/ 1// production in hadronic collisions assume that the cc pair is produced by a 2 — )■ 1 partonic process 
where both initial particles are two gluons carrying some intrinsic transverse momentum kj. The sum of the 
gluon intrinsic kj is transferred to the cc pair, thus to the 7/1// since the soft hadronisation process does not 
significantly alter the kinematics. This is supported by the picture of the Colour Evaporation Model (CEM) at 
LO (see [14] and references therein) or of the Colour-Octet (CO) mechanism at [15]. In such approaches, 
the transverse momentum pf of the // 1// comes entirely from the intrinsic kf of the initial gluons. We will 
refer to this production mechanism as to the intrinsic scheme. 

However, this is not sufficient to describe the pj spectrum of quarkonia in hadron collisions. Recent theo- 
retical works incorporating QCD corrections or i-channel cut contributions have emphasised [16, 17, 18] that 
the Colour-Singlet (CS) mediated contributions are sufficient to describe the experimental data for hadropro- 
duction of both charmonium and bottomonium systems without the need of CO contributions. For instance, 
as illustrated by Fig 1, the yield predicted by the LO CSM [19] reproduces correctly the PHENIX, CDF and 
ALICE measurements without resorting to any colour-octet mechanism nor parameter fitting. Furthermore, 
recent works [20] focusing on production at e+e~ colliders have posed stringent constraints on the size of CO 
contributions, which are the precise ones supporting a 2 — 1 hadroproduction mechanism at low p-p [14]. 

As a consequence, J /^i production at low and mid pr likely proceeds via a 2 — > 2 process, which we 
refer to as the extrinsic scheme, such ?&g + g^J/\^ + g, instead of a 2 — > 1 process. The former 2 — !■ 2 
kinematics is then the most appropriate to derive CNM effects at RHIC, and to provide predictions at LHC 
energy [21, 22]. One could also go further and consider more than two particles in the final state, as expected 
from the real-emission contributions at NLO and NNLO [16]. It is clear from the yield polarisation [23] that 
these contributions start to dominate for pr above 1 — 2mc. The effect of more partons in the final state is to 
increase the difference between the results obtained in both schemes. However the implementation of NLO 
and NNLO codes in a Glauber model with an inhomogeneous shadowing is not yet available. 

In this work, we present our results for the rapidity and transverse-momentum dependence of the nuclear 
modification factors, RdKu and /?cp obtained using the extrinsic scheme for different collision centralities. We 
compare them with the new PHENIX data [9]. 
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Fig. 1: dOjJ^^ /dy\y^Q x Br from gg fusion in pp collisions for ^/s from 200 GeV up to 14 TeV compared to 
the PHENIX [24], CDF [25] and ALICE data [26, 27]. 



2 Our approach 

In order to describe J /^i production in nuclear collisions, our Monte Carlo framework [11, 1 2] is based on the 
probabilistic Glauber model. The nucleon-nucleon inelastic cross section at y^^Aw = 200 GeV is taken to be 
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(^NN = 42 mb and the maximum nucleon density to be po = 0.17 nucleons/fm^. We also need to implement 
the partonic process for the cc production model that allows us to describe the pp data and the CNM effects. 

For a given J/y/ momentum (thus for fixed rapidity y and pr), the processes discussed above, i.e. the 
intrinsic g + g^cc^J/\lf{ +X ) and the extrinsic g + ^ —S- 7/ 1// + will proceed on the average from initial 
gluons with different Bjorken-x. Therefore, they will be affected by different shadowing corrections. 

In the intrinsic scheme, the measurement of the //i// momentum in pp collisions completely fixes the 
longitudinal momentum fraction of the initial partons: 

■^1,2 = -^exp(±y) =x?,2(3'>/'r), (1) 



with nij = y +pj,M being the //i// mass. 

In the extrinsic scheme, the knowledge of the y and pj spectra is enough to fix x\ and X2 . Actually, the pres- 
ence of a final-state gluon introduces further degrees of freedom in the kinematics, allowing several (x\,X2) 
for a given set {y,PT)- The four- momentum conservation explicitly results in a more complex expression of 
X2 as a function of {xi,y,pT): 

xiniT^/s^e-y -M^ 
X2 = , , , ^ • (2) 

Equivalently, a similar expression can be written for xi as a function of (x2,y,/?r). Models are then mandatory 
to compute the proper weighting of each kinematically allowed {xi^X2). This weight is simply the differential 
cross section at the partonic level times the gluon PDFs, i.e. g{xi , iJ.f)g{x2, IJ-f) do^^^j /y,^g/dy dpj dxidx2. 
In the present implementation of our code, we are able to use the partonic differential cross section computed 
from any theoretical approach. In this work, we shall use the Colour-Singlet Model (CSM) at LO at LHC 
energy, which was shown to be compatible (see Fig. 1) [17, 19] with the magnitude of the /^j- -integrated 
cross-section as given by the PHENIX pp data [24], the CDF pp data [25] and the recent LHC pp data at 
/y]^ = 7 TeV [27] and ^/m^ = 2.76 TeV [26]. 

To obtain the yield of 7/ 1// in pA and AA collisions, a shadowing-correction factor has to be applied to the 
J/y/ yield obtained from the simple superposition of the equivalent number of pp collisions. This shadowing 
factor can be expressed in terms of the ratios Rf of the nuclear Parton Distribution Functions (nPDF) in a 
nucleon belonging to a nucleus A to the PDF in the free nucleon: 



f^{x,Q^ 



Rf(x,Q^)= ' ^ , i = q,q,g. (3) 

I \ J J f nucleon f ' v/ 



The numerical parameterisation of Rf{x,Q^) is given for all parton flavours. Since quarkonia are essen- 
tially produced through gluon fusion at RHIC [14], we restrict our study to gluon shadowing. Several shadow- 
ing parametrisations are available. Here we will consider three of them: EKS98 [28], EPS08 [29] and nDSg 
at LO [30]. Recently, a new parametrisation with fit uncertainties, EPS09 [31], has been made available. Yet, 
in the case of gluons, the region spanned by this parametrisation is approximately bounded by both the nDS 
and EPS08 values. The central curve of EPS09 is also very close to EKS98. We consider sufficient to use only 
EKS98, EPS08 and NDSg. The spatial dependence of the shadowing has been included with a shadowing 
proportional to the local density [32, 33]. 

The second CNM effect that we take into account concerns the nuclear absorption. In the framework of 
the probabilistic Glauber model, this effect is usually parametrised by introducing an effective absorption 
cross section dabs- It reflects the break-up of correlated cc pairs due to inelastic scattering with the remaining 
nucleons from the incident cold nucleus. Here we choose four values of the effective absorption cross section 
(Cabs = 0,2.8,4.2,6 mb) following our previous works [12, 13]. 



3 Results 

3.1 Rapidity and transverse-momentum dependence of Rd Au 



We first present our results for the nuclear modification factor Rjau which characterises the J/yf suppression 
in dAa collisions. It is the ratio obtained by normahsing the J /iff yield in dAu collisions to the J /iff yield 
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in pp collisions at the same energy times the average number of binary inelastic nucleon-nucleon collisions 




(a) EKS98 (b) EPS08 (c) nDSg 



Fig. 2: J/y nuclear modification factor in dAu at ^/snn = 200 GeV for peripheral (upper part) and central 
(lower part) collisions. The four curves correspond to different values of the effective absorption cross sec- 
tion using different gluon shadowing parametrisations: (a) EKS98, (b) EPS08, (c) nDSg. The data are from 
PHENIX' [9]. 



In Fig. 2, we show RdAu vs y obtained for different shadowing parametrisations, EKS98, EPS08 and nDSg. 
We focus only on the extrinsic scheme. Our curves are compared to the PHENIX data [9]. The lower panels 
in each of Fig. 2 (a), (b) and (c) refer to central collisions (centrality class: 0-20 %, i.e. the 20 % most central 
collisions) and the upper panels to peripheral collisions (centrality class: 60-88 %). Our previous study [13], 
based on older PHENIX data [8] suggested that the effective absorption cross section which reproduced the 
most accurately the data was Gabs ~ 3 — 4 mb. Among the four different values Gats = 0,2.8,4.2 and 6 mb, 
which we have been considered here, the best match seems to be between 2.8 and 4.2 mb. The agreement 
is good for the most central collisions. For peripheral ones, none of the gluon nPDFs which we used is able 
to accommodate with the most backward data. We also note that the precision of the data does not allow 
to distinguish between the different shadowing parametrisations. These results show similar features to those 
Ref. [9], where a„/„ is taken to be 4 mb. This was expected since EPS09 shadowing is approximately bounded 
by EPS08 and nDSg and its central curve is close to EKS98. 

We now turn to the discussion of the transverse-momentum dependence of 7?^au in the mid-rapidity region. 
Once more, we would like to emphasize that it can only be predicted if one works in the extrinsic scheme. 
Our results for different centrality classes for EKS98 are shown on Fig. 3, for EPS08 on Fig. 4 and for 
nDSg on Fig. 5. RdAu is found to increase with pj. This is due to the increase of X2 for increasing pr which 
follows from Eq. 2. This effect is more pronounced for the EPS08 than for EKS98 and nDSg due to its 
stronger antishadowing. Note that the centrality dependence induced by the anti-shadowing - via its strength 



Note that the PHENIX points showed here do not include a global systematic uncertainty of ±10% for the peripheral data 
and of ±8.5% for the central ones. 
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dependence on the local nuclear density - is increasingly compensated by that of the break-up probability 
for increasing Oabs- For central collisions, the production point can be well inside the gold nucleus where the 
anti-shadowing is expected to be stronger but where the break-up probability is also larger. 

Our results are also compared to the most recent PHENIX data [34] which suffer from rather large exper- 
imental uncertainties for increasing pj. The agreement with the data is reasonable. In addition to the plot of 
^rfAuVS Pt in the mid rapidity region of PHENIX, we show in the appendix our predictions for backward and 
forward rapidities to be compared to forthcoming data. 




Extrinsic pT, EK59S shadowing 

°,bs - 4 2 mb 




Fig. 3: //)//■ nuclear modification factor in dAu at ^/snn = 200 GeV vs pf with \y\ < 0.35 in different centrality 
classes for 4 effective absorption cross sections using the EKS98 gluon nPDF. 
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Fig. 5: Idem as the Fig. 3 for nDSg. 



3.2 Rapidity dependence of Rcp 



In this section, we will discuss the rapidity dependence of /?cp which give specific information on the cen- 
trality dependence of the CNM. This quantity has the advantage to be a ratio in which many of the systematic 
uncertainties of the data cancel. It is the ratio between central and the peripheral RdAu, 
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(a) EKS98 



(b) EPS08 



(c) nDSg 



Fig. 6: Rcp nuclear modification factor in liAu collisions at ^/sfj^ = 200 GeV versus y for 4 values of Gats 
using: (a) EKS98, (b) EPS08, (c) nDSg. The data are from PHENIX" [9]. 



Fig. 6 presents our results for Rqp versus y for three gluon nPDFs (EKS98, EPS08, nDSg) and the same 
four values of Oats as above. We have already discussed the corresponding preliminary data from PHENIX 
in [13], from which we performed fits of the effective break-up cross section. We had shown at that time that 

^ Note that the PHENIX points do not include a global systematic uncertainty of ±8.2% . 
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a constant value of Oabs was acceptable when the EPS08 nPDF was chosen. As we obtained in our complete 
fit [13] taking into account all types of experimental errors [13], the comparison with the published PHENIX 
data shown on Fig. 6 suggests a Oabs smaller than what would be expected from the comparison with RdAu 
presented in the previous section. Our curve for EPS08 seems to better reproduce the most forward points, 
while it slightly misses two of the three mid-}' points. A strong shadowing seems in any case needed to account 
for these data. 



4 Conclusions 

We have evaluated the rapidity, the centrality and the transverse-momentum dependence of Cold Nuclear Mat- 
ter effects -essentially the shadowing- on J /iff production versus rapidity and transverse momentum in dAa 
collisions at Onn = 200 GeV and compared our predictions with the latest PHENIX data. We have used our 
probabilistic Glauber Monte-Carlo framework, JIN, which allows us to encode 2 — !■ 2 partonic mechanisms 
ioT J/y production. In particular, we have used the CSM at LO which is now recognised to correctly account 
for the bulk of the J/\j/ cross section in pp at RHIC. 

We have used three gluon nPDFs (EKS98, EPS08 and nDSg) and considered a reasonable range of effec- 
tive absorption cross sections, a„bs = 0,2.8,4.2,6 mb. Our results, compared to the most recent PHENIX data 
[9] are in agreement with our previous study [13] where Gabs ~ 3 — 4 mb was suggested from the compari- 
son with RdAu and Gabs ~ 2 — 3 mb from the comparison with /?cp- This difference may have some physical 
meaning but the uncertainties both in the knowledge of gluon (anti-)shadowing and in the experimental data 
preclude drawing any strong conclusions. Finally, we reassess that EPS08 with a strength proportional to 
the local nuclear density combined with a 2 — 2 partonic process is found to reproduce fairly well the most 
forward Rqp data. 
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A Appendix: RdAvy^ Pt for different rapidities and centrality classes 

In addition to the plot of RdAu^^ Pt in the mid rapidity region of PHENIX, we show in this appendix our 
prediction for backward and forward rapidities. 




2468 2453 24 



(GeV/c) (GeV/c) p^ (GeV/c) 

(a) EKS98 (b) EPS08 (c) nDSg 

Fig. 7: J/xjf nuclear modification factor in dAu at y^IjvTv = 200 GeV vs pr integrated on the centrality, for four 
effective absorption cross sections using a) EKS98, b) EPS08, c) nDSg in the 3 rapidity regions covered by 
PHENIX 
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(a) EKS98 (b) EPS08 (c) nDSg 

Fig. 8: Idem as the Fig. 7 for the centrality class 0-20 %. 
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Fig. 9: Idem as the Fig. 7 for the centrality class 20-40 %. 



10 



-"^o.'^<>.-o-H>-o-o<ho-o- -o-o_f 



1 — 



— 0.5 - 



P, (GeV/c) 



(GeV/c) 

(a) EKS98 (b) EPS08 

Fig. 10: Idem as the Fig. 7 for the centrality class 40-60 %. 
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